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Abstract. A quantum theory is developed for a difference-difference system which 
can serve as a toy-model of the quantum Korteveg-de-Vries equation. 



■-■ Introduction 
in ■ 
G> ! 

Q"^ , This Letter presents an example of a completely integrable 'discrete-space-time quantum model' whose 
' Heisenberg equations of motion have the form 

\D] =A0(r,n)(^(T-l,n)+ (/)(T-l,n)0(r-l,n-l). (1) 

T-H . By 'discrete... model' we mean 
>■ 

! (i) an algebra 'of observables' 

I whose generators 0„ are labeled by integer numbers n which are regarded as a (discrete) spatial 

' variable; together with 

in : 

On ■ (ii) an automorphism Q, 

'. 

. whose sequential action 
CIh' X (0) = X e $ 

D ; 

^ ■ Q : . . . 1-^ X (r — 1) 1-^ X (r) 1-^ X (r -I- 1) 1-^ . . . 

• 5^ I is viewed as the (discrete) time evolution. Thus, we intend to produce a pair <I>&Q such that the 

^ ■ evolution <pn{T) of generators, in the natural notation 

S ■ 

■ - - ' 4){T,n)= 0„(r), 

obeys the system (1). 

Complete integrability is understood as the existence of a commutative subalgebra 'of conservation 
laws' preserved under time evolution and spanning, in a sense, half of the algebra of observables: it 
is commonly believed that a Hamiltonian system may be either 'completely' nonintegrable possessing 
only a few conservation laws due to its manifest symmetries, or completely integrable enjoying a whole 
lot of conservation laws, one per degree of freedom. The commutative subalgebra which we encounter 
in this Letter definitely contains a lot of conservation laws but the question of how many is left to be 
answered elsewhere. 

Actually, we deal here not with a single model but rather with a family of them <I>&Q(A), each model 
being related to a certain value of a complex parameter A in (1). Moreover, all their evolution automor- 
phisms Q(A) turn out to be mutually commuting and sharing the common subalgebra of conservation 
laws. 

The order of presentation is as follows. In Section 1 we introduce an algebra of observables which is 
basically the same lattice U{1) exchange algebra which appeared already in [FV93, 94]. Naturally, the 
behaviour of that algebra depends on the value of a constant q involved in the commutation relations. 
For simplicity we shall assume that q is a root of unity. 



*On leave of absence from Saint Petersburg Branch of the Steklov Mathematical Institute, 
Fontanka 27, Saint Petersburg 191011, Russia 



After some preliminaries (Sections 2 and 3) we pick in Section 4 from the algebra of observables a 
set of 'Fateev-Zamolodchikov i?-matrices' Rn(A) which satisfy a chain of Yang-Baxter equations 

R„_l(A) Rn(AM) Rn-l(/^) = Rn(/^) Rn-l(A/x) Rn(A) 



amounting to mutual commutativity 



Q(A)Q(m)= Q(/x)Q(A) 



of their properly defined 'ordered products'. The family Q (A) provides the demanded commuting 
(inner) evolution automorphisms 

(j){T,n) = Q-^(/)„Q^ 

and doubles as their common conservation laws. 

In Section 6 we eventually establish that these evolutions do solve the equations (1). Prior to that, 
in Section 5, we discover one more face of the family Q (A). As a function of A it proves to satisfy the 
Baxter equation [Bax] 

q^^' Q (A) t (A) = a^(A) Q (q'U) + (5^(A) Q (qA) 

which in turn makes Bethe ansatz equations to emerge in a purely algebraic context. 

While the quantum system (1) seems to be a recent invention (it has been looked at, albeit from a 
somewhat different angle, in [FV92, 94]), its classical counterpart has been around for quite a while. It 
was introduced (in a somewhat different form) by Hirota back in 1977 [H] as an integrable difference- 
difference approximation of the sine- Gordon equation but eventually proved far more universal making 
perfect sense as a lattice counterpart of numerous integrable equations including that of Korteveg and 
de Vries. To conclude the Introduction we shall list various continuous limits and alternative forms of 
the (classical) system (1). This should give some idea of where our model fits into the scheme of things 
accepted in Soliton Theory. For more of that and a comprehensive list of relevant references see [NC] . 

Pairs of integers (r, n) may naturally be viewed as vertices of a plane lattice. In the shorthand 
notation 




for a quartet of vertices enclosing some elementary cell of that lattice, with subscripts instead of paren- 
thesized arguments and without bold letters reserved for the quantum case, the classical equations (1) 
read 

(f^B^A - ^^^C - A (0B^D - 0A^c) = 0- 

• At least one continuous limit is already quite apparent. Let us put the lattice on the coordinate 
plane (t,x) in such a way that vertices (r, n) go to points (A~^Ar, An): 



A 



B 



D 



t 

A/A 



If now one manages to find a family of solutions to (1) depending on the lattice spacing A and tending to 
a smooth function v{t, x) as A goes to zero then that function can be easily seen to satisfy the equation 



This equation hardly needs any comment although at this point the use of a nonlinear difference- 
difference equation to model a linear differential one seems difficult to justify. Once, however, one takes 
it as a Minkowsky version of the Cauchy-Riemann equations things start to look like a unified approach 
to conformal invariance and integrability. 

• Let us now perform a more sophisticated superimposition (r, n) — > ^^"^ + An): 




This time we come to the close relative of the KdV equation 



for any solution of which the potential 



h + 9xxx - 3 — - — = 



u = 



solves the KdV equation itself 



Ut + Uxxx — 6 UUx = 0. 



• To demonstrate a somewhat different scenario let us recall how to turn (1) into the sine-Gordon 
equation. Before performing a continuous limit we switch in (1) to the function <f{T,n) 

(/.(r,n) = e(-i)"*^(^'") 

that brings (1) back to its original form 

sini((^^ -ip^-ifc +<^^)+Xsm^{(p^+ip^ + ip^ + ip^) = 0. 
Let the lattice cells now look like 



A 



C 



B 



D 



A 



^a4 



and rescale the constant in the equation making it dependent on the lattice spacing 

If now for a fixed value of the constant m one finds a family of solutions tending, as A ^ 0, to a 
smooth function (p{^,rf) then this function satisfies the sine- Gordon equation 

2 

Of course the original function (j) can not survive under this continuous limit becoming badly oscillating. 

• Although the mere ability of (1) to unify KdV and SG equations makes it a reasonable prospect, 
one is still left to wonder why the would-be universal lattice equation does not look special enough for 



that. And indeed there exists a far more spectacular version of (1). One can try to guess its form 
recollecting that the KdV equation looks best in its Krichever-Novikov reincarnation^ 



ft + S[f]U = 
where S[f] stands for the Schwarz derivative 

fxxx 3 / fxx ^ 



S[f] 



fx 2 \ /a; 



Any decent lattice approximation for the /-equation should employ the cross-ratio as a difference coun- 
terpart of the Schwarz derivative and it is not difficult to spot the right one: 

if A ~ /b)(/c ~ /d) ^ const. 

if A ~ /c)(/b ~ /d) 

Indeed, this equation does produce, provided const= A~^, in the above 'parallelogram' continuous limit 
exactly the /-equation being as well a completely integrable difference-difference model in its own right. 
And just as the continuous /- and (/>-equations are tied up by the Fuchs formula 

their lattice counterparts are connected by the map defined by 

~r{fn ~ fn~l) = — — 
A (Pn-l(Pn 

where we omitted the argument r common for all entries and moved the remaining one to the subscript 
position. 

The cross-ratio version will find extensive use in the forthcoming paper addressing higher-order 
equations (KdV hierarchy), equations with two fields (NLS hierarchy) and 2-|-l-dimensional equations 
(KP hierarchy). Unfortunately, despite of its virtues this version has not yet been really useful in the 
quantum theory where at present we are only able to handle a free-field sort of algebra of observables 
associated with the (/)-equation. 

• By the way, at this point one finds himself well prepared to design a lattice version of the u- 
equation. Introducing 

(Pn-l + 9n+l 

SO that 

^ = 4 (/"+! ~ fn){fn-l - /n-2) 

(/n+1 — fn-l){fn — fn-l) 

and 

one eventually transforms (1) into 

K - '"d)K - ^c) + - 1)K^B - «c^d)^ = 

approximating the KdV equation in its original form. 

^As this is already the third form of the KdV equation we have met, so far, and the arrival of a fourth one, the so-called 
modified KdV equation, is imminent, let us recall how all these forms interact 

/ ft + s[f]u = 

I 

- ^ + <t)... - 3%^ = 



2 fa: 



Pt + Pxxx -Gp Px =0 



U 



P^ + Px = = —\S[f\ Ut + Uxxx ~ 6uUx ^ 



1 Algebra of observables 

All the way through this Letter a constant q will be an odd root of unity 

q2^+l = 1 

and an integer number A^(> 3) called the spatial period will be odd. The only algebra of observables 
in use will be the 'exchange' algebra $ with generators (f)n subject to 



(i) commutation relations 
with 



complemented by conditions 



(iii) 



e(n) = 1 n = 1,3,5, ...,iV-2 
e(n) = n = 0,2,4, ...,iV-l 
e{n + N) = e(n) + 1, 



0r^ = l, 



(f) (f) n+]\f does not depend on n. 



The simplifying condition (ii) is natural if not really essential. In (iii) one easily recognizes a quasiperi- 
odic boundary condition leaving in the algebra of observables only independent generators0, for 

instance, (fyo, 

That algebra contains a useful 'current' subalgebra W with generators 

0n+l 
W„ = —r 

that are easily seen to obey 
(i) commutation relations 

w„,_i w„ = q^ w„_i 
^ n = ^ m \m — n\ ^ 1 (mod A^), 

and conditions 
(ii) 

(iii) 

Of course, (iii) is just the periodic boundary condition. The whole algebra $ is, in a sense, one degree of 
freedom larger than that current subalgebra which has A'^ independent generators and a central element 

C = q( 0n+Ar)^ = qWiW3... WArW2W4... WAr_i. 

A remark is in order here. The algebra of observables was designed with the 'second' periodic KdV 
bracket (aka the Virasoro algebra) 

i{n(x), u{y)} = 2iu{x) + u{y))6'{x - y) - 6"'{x - y) 



*Af+l is even, all right. 



in mind. Indeed, that bracket corresponds to the sign-function-bracket (the exchange algebra)^ 

(l){y)} = -\ s\gD.{x - y)4,{x)4,{y) 
in the i?!)-language and the delta-prime-function-bracket (the current algebra) 

^{p{x),p{y)} = 5'{x - y) 

in the p-Ianguage (see Introduction) . Both are easily seen to be classical (q = e*^''' , ^ ^ 0) continuous 
{4>iat ~ 4'cont-, wiat ~ 1 + 2Apcont) limits of the above commutation relations. A natural question 
arises whether the Virasoro algebra itself has a reasonable lattice counterpart. The answer seems to be 
affirmative [FT, V, B, Fe] but this is another story. 

2 Fateev-Zamolodchikov i?-matrix 

Let two operators u and v satisfy the condition 

U2'^+1= V 2^+1 = 1 

and obey Weyl's commutation relation 

u V = q^ V u . 

Then, as was found in [FZ], the pair of functions of complex variable taking values in the algebra 
generated by u and v 

i?i(A) = r(A, u) i?2(A)= r(A, v), 

where 

i 

k=-t 

k I 
Pfc(A) = p_fc(A) = q^'' n(l - Aq-20-1)) n (1 - Aq2^) ^<k<l 
j=l j=k+l 

and n? • • • = n'+i • • • = 1, 

satisfies the 'braid' Yang-Baxter equation 

RiiX) R2iXfi) Riifi) = i?2(Ai) RiM i?2(A) . 

Before going on let us compile a list of some useful properties of the function r(A, z). We shall often 
use them in remaining sections, sometimes not mentioning it explicitly. In what follows the second 
argument of r is always assumed to satisfy the condition 

= 1. 

(i) r(A, z) is, up to a constant factor, the only polynomial in A of degree i satisfying the functional 
equation 

{X + z)r{X,qz) = {1 + Xz)r{X,q'^z). 

(ii) r{X,z) satisfies the functional equation 

(1 + A)(l + Aq) r{qX,z) = {1 + Xz){l + Xz'^) r{q-^ X, z). 

(hi) 

r{X,z)= r{X,z-') 



^Strictly speaking, 5{-) means here the 27r-periodic 5-function while sign(-) is a function coinciding with the usual sign 
function in the interval [— tt, tt] and extending quasiperiodically elsewhere: sign(a:: + 2n) — sign(a;) + 2. 



(iv) the product of r{X,z) and r(A does not depend on z: 

e 

r(A, z) r{X-\ z) = q{X) = const • A"^ 11(1 + Aq'^). 



(v) At the point A = the function r{X,z) turns into a 'truncated' theta-function 

e 

r{o,z) = eiz)= J2 q' 

satisfying the functional equation 



■fz" 
k=-e 



zeiqz) = eiq'^z). 

The operators 

01 = e{u) 92 = e{y) 
obey Artin's commutation relation 

616261 = 626162 
thus providing a 'free- field' model of the braid group i?2- 

(vi) 

r(l, z) = const 

Some of these statements are quite transparent, some are less so. We will present their proofs in a more 
detailed paper. 

3 Cyclic product 

Let be a free algebra^] with generators Denote 

(i) by <j the homomorphism projecting TV to 

?(x y) = c^(x) <;-(y), 

(ii) by +Ar the shift-by-period automorphism of the algebra W 

( Wn)+N = Wn+AT 
(x y)+N = X+Ary+AT. 
Assign to monomials of W 'overlap' numbers 

+00 L 
p=— 00 i,j=l 

and utilize them in the coboundary which defines another multiplication ★ in W: for basis elements it 
rCctds 

f =q2(Kfg)-Kf)-Kg))f g 
and extends bilinear ly elsewhere. Then for any two elements x , y of 

?(x*y) = ?(y* x+at). 

We omit the proof for it is a straightforward computation. 

The purpose of this <; (*) construction must be clear for those familiar with the Quantum Inverse 
Scattering Method: once we decide to do without an 'auxiliary space' and go for a purely algebraic 
version of the i?-matrix approach we need some direct method of computing what used to be auxiliary 
space traces. 

^It is essential that the algebra W, as opposed to the algebra W of Section 1, is devoid of the periodic boundary 



4 Conservation laws 



We are going to prove that 'cyclic ordered products' 

\J{X) = <; (Ri(A)* R2(A)*...* Rjv(A)) 

e 

= E ^^'"iPjW^) R2(A)R3(A)... R;v-i(A)(pfc(A)w^) 

j,k=-e 

of FZ i?-inatrices 

R„(A) = r(A, w„) Rn(A) = r(A, w„) 

commute with each other: 

U(A)U(m)= U(m)U(A). 
Indeed, the commutation relations of w 's translate into N copies of the Yang-Baxter equation 

R„_l(A) R„(A/X) Rn-l{fl) = Rn(/^) Rn-l(A/x) Rn(A) 

Rn+Af(A) = Rn(A) 

that is believed to ensure the commutativity of ordered products of those i?-matrices. It is not however 
evident how to make this idea work in the periodic case where a naive ordered product, say, R i R 2 • . . R at, 
makes little sense. The <; (*) 'product' has better chance to deliver, since, as we know from the preceding 
section, it does not at least depend on the starting point: 

... = <j (Ro* Ri*...* R N-i) = <; (^Ri-k R2* ■ ■ ■-'^ ^n) = ■ ■ ■ ■ 

And once we get the ordered product right the commutativity check becomes a matter of familiar 
i?-matrix machinery complemented by the rules of (★) 'multiplication' (and also by the definition 

r'^^(A) = R„(A-i)/^(A) ensuring that ^ ( r'^^(A) Rn(A)) = 1): 

U(A)U(/x) =^ (Ri(A)* R2(A)*...* R;v(A)) ^ (Ro(m)* Ri(m)*...* Riv-i(Ai)) 
= ? (Ri(A)* Ro(/i)* R2(A)* Ri(/i)*...* Riv(A)* Rn-i{i^)) 
= ^ (Ro'(^)* Ro(^)* Ri(A)* Ro(/^)* R2(A)* Ri(/x)*...* RAr(A)* Riv-i(^)) 
= ^ (Ro(^)* Ri(A)* Ro(Ai)* R2(A)* Ri(/x)*...* Rjv(A)* Riv-i(/x)* R]^^(^)j 
= ? (^Ri(/i)* Ro(A)* Ri(^)* R2(A)* Riifi)*...^ RnW* Rjv-i(/u)* Rn\^)^ 
= ? (Ri(/x)* Ro(A)* R2(/x)* Ri(A)* R2(^)*...* Rjv(A)* Rjv-i(/x)* 



= ^ (Ri(/i)* Ro(A)* R2(/i)* Ri(A)*...* Riv-i(^)* Riv(A)* Riv-i(/x)* rV(J)) 

= q (Ri(/i)* Ro(A)* R2(m)* Ri(A)*...* Rjv(Ai)* Rjv-i(A)* Riv(^)* rV(^)) 
= q (Ri(/x)* Ro(A)* R2(/x)* Ri(A)*...* Rjv(/x)* Riv-i(A)) 
= ^ (Ri(m)* R2(m)*---* Riv(/x)) ? (Ro(A)* Ri(A)*...* Riv-i(A)) = U (/x) U (A). 
We conclude the Section with two remarks. First, to make formulas easier on the eyes we shall adopt 

o 

the notation ... for 'cyclic ordered products' like the ones above, for instance, 

o 

U(A)=n ^niX) 

o 



Second, for a reason that will become apparent in Section 6 the family U (A) should be properly 
'normalized' 

Q(X)- 

that, of course, does not spoil its commutativity 

Q (A) Q (i^l) = Q (^) Q (A) 

and polynomiality in A. 

5 Baxter equation 

In the Quantum Inverse Scattering Method language the family Q (A) would be called the fundamental 
transfer-matrix [TTF] as opposed to the usual nonfundamental one which on this occasion has the form 
[G, V] 

o 

t(A) = trn(W„L(A)) 

where the matrices W^i and L{X) are 

W„ = I ™" _l I with w| = w„ 

tr denotes the matrix trace and * (hidden in the product symbol) combines what it used to be with the 
standard matrix multiplication. In the decyphered form it reads 

t (A) = A'^ ^ q2A:ifc]v ^2(A:iA:2+fc2A:a+...+fcjv-ifc]v+A:jvfci) ^fci ^te 

fcl,fc2,...,fcjv=±5 

with 

q2 = q ^ 

In other words, we have another polynomial in A, this time of degree , which, as we know from past 
experience, 

(i) commutes with itself 

t(A)t(/.)= t{^,)t{X), 

(ii) commutes with Q (A) 

t(A)Q(^)= Q(^)t(A) 

and, as we are going to see, 

(iii) satisfies together with Q (A) the Baxter equation [Bax] 

q^^' Q (A) t (A) = a^(A) Q (q-^A) + 5^(A) Q (qA) 

with 

a(A) = A - 1 
5(A) =qA + l. 

The first two items of this list are actually superseded by the much stronger third one the proof of which 
may go as foUowsFl: 



"To be precise, only odd degrees are present. 
**The first proof of (iii) for the model in question was obtained by R. Kashaev [K] while the very idea that the fundamental 
transfer-matrix Q (A) can serve as a Baxter's Q-operator probably belongs to E. Sklyanin [S]. See also [PG] addressing 



o o 

U (A) t (A) = tr n ( Rn(A) W,_iL(A)) = tr [] C„(A) 



with 



C„(A) = B-^ W;ii R„(A) W„_iL(A) \N^B, 
B -- 



1 1 
1 



Using the obvious relation 

one gets 
with 
C{X,z) 



Wn-i r{X, w„) w„^i = r(A,q w„) 
C„(A) = C(A, w„) 

^ z^{Xr{X,q~^z) — r{X,qz)) {Xz^ + z~^) r{X,q~^z) — {z^ + Xz~^) r{X,qz) \ 



V 



22 r{X,qz) 



{z-2 + Xz 2) r{X,qz) 



I 



The first property of the function r(A, z) (see Section 2) says that the upper off-diagonal 
element of this matrix vanishes (provided z^^^^ = 1 and z^ = z"^). This yields immediately 

o o 

U (A) t (A) = n w„) + n ^^(A, w„) 



where a and d denote the diagonal elements of the matrix C: 

C{X,z) 

It remains to verify that 



a{X,z) 
* d{X,z) 



aiX,z)=q-' (A-l)r(q-iA,z) 

d{X,z) =q"^'(qA+l)r(qA,z). 

We omit this part of the proof for it is neither difficult nor instructive. 

To conclude, let us recall what use might be made of the Baxter equation. The function Q (A) is a 
polynomial in A with coefficients coming from the subalgebra of conservation laws. This allows, probably 
at the expense of considering a proper completion of the algebra of observables, the introduction of 
commuting 'roots' of Q (A) 

Q(A).n(i-A) 

""j rjfc = rfe r j. 
The Baxter equation says, in particular, that 

a^(A)Q(q-iA) + 5^(A)Q(qA) =0. 
This is just the famous system of the Bethe ansatz equations which look more familiar in the form 



Sirk) 



N 



6 Equations of motion 

We are going to see that the generators of the algebra of observables evolve 

0(A|r,n)= Q-^(A)(/>„Q^(A) 

according to the equations (1) which in ultimately accurate form read 

(A I r, ra) (A I r, — 1) + A (A I T, n — 1) (A I T — 1, n — 1) 

= A (/> (A I T, n) </) (A I T — 1, n) + (/> (A I T — 1, n) (A I r — 1, n — 1). 
We shall cover the distance in four short steps. 

(i) Commutation relations between (f) 's and w 's have the form 

(t>m^n= '^n4>m m ^ n {uiod N) . 

This prompts us to read the first property of the function r from Section 2 as 

{ (p n+l (p n + >^ n (f> n-l) Rn(A)= R „(A) (A „+i n + (f) n (p n-l) ■ 

(ii) Viewing this equation as a fragment of the 'big' one 

e 

^^^'^ (Pj^m+l) Rm+2-- - Rn-l(0n+10n + A0n</>n-l) RnRn+1--- Rm+AT-l {Pk^L+N) 

j,k=-e 

e 

= 'i^^'' (Pj^ii+l) ■ ■ ■ Rn-lRn (A </>n+l n + 4> n 4> n-l) Rn+1 • • • Rm+Af-1 (pk^m+N) 

j,k=-e 

and pulling </> 's to the outside(s) we get 

U (A) + A (/)„ U (A) = A (?!)„+i U (A) + U (A) (pn (pn-1- 

(iii) At the point A = it turns into 

</)nU(0)= \J{0)(pn<Pn-l 

suggesting more subtle 

(pn^iO) = U(0)(/)„„i. 

As a matter of fact, the latter does hold. We omit the proof for it is too case-specific. 

(iv) (ii) and (iii) combined yield 

</>n 0n-l Q (A) + A Q (A) (/>n-l = A (/>n Q (A) 0„ + Q (A) (p n (p n-l 

which is nothing but (1) with cut away common factors Q ""^(A) . . . Q'^~^(A). 

So, we have met the last objective of the Letter. We conclude it with two remarks. 

• A consistent approach to the subject should probably distinguish between observables and their 
automorphisms rather than mix them up as we did in this Letter. It would be only natural to deal not 
with the i?-matrices but directly with automorphisms they represent: 

TZn{X) : A+q w„ 

m ^ n (mod N) 
Indeed, it would follow straight from the above definition that 



and 

Un-i{x) oTZn{Xn) o7^„_l(/x) = 7^„(/i) o7^„_l(A;u) o7e„(A). 

These two relations are 'weaker' (but just sufficient!) substitutes for the two cornerstones of the whole 
scheme, which are item (i) of this Section and the Yang-Baxter equations of Section 4. Why then care 
whether those 'i?-automorphisms' are inner or not? They happen to be inner in our particular case 
but even there some other important automorphisms are outer anyway [FV93]. In other cases one pays 
a dear price for 'inclusion' of i?-matrices in the algebra of observables [BR, F, FV95]. Unfortunately, 
some difficulties of the 'automorphism' approach made the author to choose for this Letter the more 
familiar 'inner' route. 

• It would be useful to know whether the equations of motion (1) provide the exhaustive information 
about the model. In other words, is <p (r, n) the only solution to the Cauchy problem 

4>{T,n)\r=0 = 4>n 

for the system (1)? Not quite, and it is easy to see why. First, in our solution the 'quasimomentum' c 
(see Section 1) does not evolve but the equations (1) do not know about that. Also, they do not feel 
whether we multiply R„(A)'s 'from left to right' or R„(— A)'s in the opposite order. As a matter of fact, 
there are no other sources of nonuniqueness. So, the local structure of (1) is good and all one needs to 
achieve the ultimate uniqueness is a couple of extra global conditions. This point will be described in 
more detail elsewhere. 
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